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Superradiance, subradiance, and suppressed superradiance 
of dipoles near a metal interface 
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We theoretically characterize the collective radiative behaviour of N classical emitters near an 
interface between different dielectrics that supports the transfer of surface plasmon modes into the 
far-field of electromagnetic radiation. The phenomena of superradiance and surface plasmons can be 
combined to amplify the emitted radiation intensity S as S = AN^So compared to a single emitter's 
intensity So in free space. For a practical case study within the paper A = 240, compared to ^ = 1 
in free space. We furthermore demonstrate that there are collective modes for which the intensity of 
the emitted radiation is suppressed by two orders of magnitude despite their supperadiant emission 
characteristics. A method to control the emission characteristics of the system and to switch from 
super- to sub-radiant behaviour with a suitably detuned external driving field is devised. 
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INTRODUCTION 



Observing light spontaneously emitted from a system 
of atoms or molecules is a versatile tool to determine the 
system's dynamics and has a number of important ap- 
plications. In atomic physics it can generally be used 
to learn about the state of atoms and molecules [l] . In 
quantum information it helps to detect entanglement be- 
tween ions or atoms [2, 3 . In biophysics it is used to iden- 
tify and track chemical compounds in biological tissues. 
While spontaneous emission is well understood for iso- 
lated emitters, emission from a system of emitters varies 
extensively with their interaction, depending on the state 
of nearby emitters and the surrounding media. Only by 
understanding the role of interactions in light emission 
can the state reliably be determined. The process of the 
collective decay rate being coherently increased is known 
as superradiance [4 , and has received renewed interest 
with the development of nanotechnologies [5 . In quan- 
tum optics, superradiance can be used to create entangle- 
ment between atoms and light [6]. Entanglement within 
a superradiant system is also being used to understand 
phase transitions in quantum systems [7 . For this rea- 
son there is strong interest in methods that affect the 
interactions in cooperative emission. 

A simple method to influencing the interaction 
strength is to place the system near a metal interface 
that supports surface plasmon modes [8^. Surface plas- 
mons are charge density waves confined to a very small 
region close to the interface. The field of surface plas- 
mons are consequently very large. This effect has long 
been used in techniques such as surface enhanced Raman 
spectroscopy [9l |10]. In quantum information science, 
surface plasmon modes on a nanowire are proposed to 
create single-photon sources and transistors [TTl [12] and 
enhance the emission properties of light emitting diodes 
[13]. Furthermore, the coupling of emitters to surface 
plasmon modes allows the collective excitation of super- 
radiant surface plasmons [14 . 



In this paper we show how superradiance and surface 
plasmons can be combined to increase the interaction 
strength between matter and light by several orders of 
magnitude. We consider emitters near a planar metal in- 
terface in the Kretschmann configuration for attenuated 
total reflection [15 , see Fig. [l] In this situation, surface 
plasmons induce a greater coupling of radiation into the 
far-field for a certain direction (the surface plamon res- 
onance angle) where the flux of radiation is two orders 
of magnitude larger [16 . Furthermore, the increased ra- 
diation energy density in the near field of surface plam- 
ons leads to a stronger coupling between different emit- 
ters [17 so that superradiance is further increased. We 
characterize the collective radiative eigenmodes of a set 
of oscillating emitters and demonstrate how their emis- 
sion can be manipulated from super- to a sub-radiant 
behaviour. 

After a brief review of superradiance and surface plas- 
mons in Sec, [nj and our model assumptions in Sec.|III|we 



will consider harmonically oscillating emitters that are 
driven by a light field in Sec. IV We demonstrate that 



their collective dynamics can be changed from super- to 
sub-radiant behaviour by changing the frequency of the 
field. In Sec.|V|we study the collective decay of initially 
excited harmonic oscillators and give a detailed descrip- 
tion of super- and sub-radiant decay modes, their decay 
rates, emission pattern, and discuss how well superradi- 
ant emission is realized in these modes. Three appendices 
contain details of our derivations. 



II. BACKGROUND 

Superradiance is a cooperative effect in which the radi- 
ation emitted by N scatterers is substantially enhanced 
compared to isolated scatterers [4 . For this effect to oc- 
cur the distance between the emitters needs to be much 
smaller than the wavelength of the radiation. The field 
amplitudes of the emitters then interfere constructively 




FIG. 1: Kretschmann configuration of a tiiin metal film in- 
terface between a prism and vacuum used to couple emit- 
ters to surface plasmon modes. Region 1 consists of a prism 
(v^ = 1-51), region 2 corresponds to a thin metal film 
(v^ — 0-08 + 4.8i), and region 3 is the vacuum. 



so that the radiation intensity is increased by a factor 
of N'^ compared to a single emitter. Furthermore, light 
emitted by one of the scatterers can induce stimulated 
emission for another scatterer so that the collective de- 
cay rate is enhanced by a factor of N if all emitters are 
initially excited. These two defining properties of super- 
radiant emission are consequences of phasing and radi- 
ation reaction which can occur both in classical and in 
quantum systems fl8 . Radiation reaction is generally af- 
fected by any surrounding dielectric material. In partic- 
ular, when atom-sized dipolar emitters are placed within 
a wavelength of a metal interface, surface plasmons can 
be generated. 

Surface plasmons consist of electromagnetic fields that 
accompany electron density oscillations (the plasmon) 
on the surface of a metal. Because the surface plas- 
mon is confined within a nm-sized region on the interface 
the electromagnetic field is concentrated. This concen- 
tration depends strongly on the geometry and the di- 
electric properties of the interface. Because of the in- 
creased intensity of surface plasmon fields, an emitter 
reacts strongly with surface plasmon modes when placed 
near a metal film and the radiative properties of the emit- 
ter are significantly affected [TF. In the following we will 
study in detail how surface plasmons affect superradiance 
of emitters near a planar interface. 



III. OSCILLATING EMITTERS NEAR A 
METAL INTERFACE 

We consider the behavior of N atomic dipoles that 
are located near a thin metallic interface that supports 
surface plasmon modes, see Fig. [l] To simplify the con- 
siderations we study an infinite planar interface so that 
the prism, which is needed in experiments for phase- 
matching, can be replaced by an infinite half-space filled 
with a non-absorbing dielectric. The dipoles are located 
at positions r^^\ ...,r^^^ in region 3 above the interface. 

To describe the system of emitters we model each 
dipole as a charged classical harmonic oscillator with 
charge q^ mass m, and resonance frequency of ujq. The 
restoring force associated with the harmonic motion orig- 
inates from the attraction by a stationary charge distri- 



bution of opposite sign (the "nucleus") located at the 
point r^'^\ Charge i oscillates about the point r^'^^ with 
an amplitude x^^^ that is of the order of Bohr's radius 
and much smaller than the individual separation distance 
between emitters Ax. The equation of motion for the 
harmonic oscillators affected by an external electric field 
E{r^t) is derived in Appendix [A| 



N 



;(i) 






m 



E{r^^\t)^^E^\r^^\t) . (1) 



The term in parentheses can be interpreted as the local 
electric field at the position r^^^ of the jth oscillator. It 
consists of the external applied field E{r^^\t) in region 
3 (e.g., a laser beam) plus the superposition of the fields 
E^'^\r^^\t) that are created by emitter i. Due to the 
presence of the interface the latter part of the contribu- 
tion is more involved and is conveniently expressed by 

/CX) 
dt'G{r,r^'\t-t')x(i\t'), (2) 
-CX) 

The dyadic retarded Green's function G describes the 
propagation of electromagnetic fields in the presence of 
the interface and takes into account the corresponding 
boundary conditions. It is frequently used to describe 
radiative systems with boundary conditions [2QH22] de- 
spite that fact that its explicit form is rather involved 
(see Appendix IbI). This is because the Green's function 
contains all information about the propagation of elec- 
tromagnetic waves between two arbitrary points in the 
presence of the interface, which is necessary to analyze 
emitters at different positions. Alternatively one could 
use an expansion of the electromagnetic field in terms of 
radiative field modes [23 which has been shown to be 
equivalent [24 . 

Within the Green's function formalism, the effect of 
surface plasmons are included through the Fresnel coeffi- 
cients, i.e., the complex ratio of the reflected (i^™i) or 
transmitted (T^™i) electric field and the incident field, 
where i is region in consideration and i — 1 the neigh- 
bouring region. Here "TM" refers to transverse mag- 
netic polarization which is necessary to generate surface 
plasmons. The Fresnel coefficients appear explicitly in 
the Green's function (see Appendix [B|); surface plasmons 
generate characteristic resonances in these coefficients. 

The dispersion relation for a surface plasmon at a sin- 
gle vacuum- metal interface. 
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can be derived from the condition that there is no re- 
flected field, i?™ = [19^. Here /cgp is the surface plas- 
mon wave vector component parallel to the interface. For 
a single interface between only two different dielectric me- 
dia the dispersion relation corresponds to the poles of the 
Green's function in frequency space. Even though the in- 
clusion of a second interface changes the dispersion rela- 
tion, Eq. ([3| is often presented within the context of both 




FIG. 2: The amplitude response in z-direction of two emitters 
in (a) free space (b) and placed at a distance of A2; = lOOnm 
from the metal interface and separated a distance Ax from 
each other. Both emitters are driven in the z-direction by an 
electric field that is detuned by 5uj from a resonance frequency 
ofcjo = 2.36x 10^^-^ 
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FIG. 3: The z-amplitude xo,3 (a) and x-amplitude xo,i (b) 
of two emitters near a metal film (solid line) and a perfect 
mirror (dashed line). The emitters are driven by a field that 
is detuned by 5uj. They are separated by Ax = 80nm from 
each other and Az = lOOnm from the metal interface. 



single and double interfaces [25 . This approximation 
suffices to determine the surface plasmon resonance an- 
gle in attenuated total reflection for thin films. However, 
it omits complex valued contributions to the wavevector 
component /Csp[19l [26], most notably complex contribu- 
tions that would correspond to radiative coupling to the 
prism (i.e., the third medium needed in attenuated to- 
tal reflection). For the Kretschmann configuration the 
narrow resonance angle of surface plasmons results in an 
increased radiation density which can be used to increase 
the coupling between emitters and radiation by several 
orders of magnitude. It is the combination of this effect 
with super-radiance that is at the centre of our attention. 



IV. SUPER- AND SUB-RADIANT EMISSION 
BY DRIVEN OSCILLATORS 

If the emitters are driven by an external electric field 



of the form E{r,t) = E^e 



ihr ^ 



-iujt 



they will eventually 



settle into a state where each dipole oscillates with the 
frequency uj of the driving field, so that x{t) = Xoe~*^^ 
[33 . It is shown in Appendix [A| that this ansatz trans- 
forms Eq. (IT]) into the algebraic equation 04] 



Eq -e 



{u^ - uj^)dijt - lmG{r^'\r^'\uj) 



X Si 



2 2 ^ 



27T 



27r 



^G{r^'\r^^\uj) 



j^i 



X 



' 



(4) 



and forms a set of 3A^ coupled equations for the Cartesian 
components of N oscillators. The indices z, j run from 1 
to N and 1 denotes the unit matrix in three dimensions. 
To find the total electric field emitted by all N sources we 
need to solve these equations for cc*^*^ and then superpose 
the fields produced by all sources. 

Figure [2|a) shows the amplitude response xo,3, i.e., 
the 2;-component of the amplitude vector xq of the oscil- 



lators, for two nearby emitters in free space driven by an 
electric field polarized in the z direction. For separation 
distances Ax less than 50 nm the resonance frequency be- 
gins to shift considerably due to dynamic dipole-dipole 
coupling between the emitters. Because in free space the 
dyadic Green's function is diagonal with respect to ra- 
diation polarization, and because the emitters are only 
driven along the z-axis, there is only coupling between 
the z-components of oscillations. 

Placing the two emitters close to the metal interface 
leads to a dramatic change in the emitter dynamics. In 
Fig.^b) the 2:-amplitude of two scatterers that are a dis- 
tance Az = lOOnm away from the interface is displayed. 
It has a similar overall shape as the amplitude response in 
free space, but an additional narrow resonance peak ap- 
pears on the tail of the primary resonance. This narrow 
feature is more prominent in Fig. Isja) which displays the 
same quantity for a fixed distance Ax = 80nm between 
the scatterers. 

The origin of the secondary resonance is the reflection 
of light emitted by one oscillator from the interface and 
its subsequent absorption by the other oscillator. In free 
space light emitted by one oscillator necessarily has the 
same polarization as the external driving field. There- 
fore, the emitters are collectively oscillating along the 
z-direction. However, light that is reflected by the inter- 
face can have a different polarization. In the case under 
consideration it induces a coupling between the x- and 
z-components of the oscillators; mathematically this is 
related to the G13 terms in Eq. Q. 

Figure IsFb) indicates that at the secondary resonance 
the dipole orientation changes from normal to almost par- 
allel to the interface. Because the two dipoles are nearly 
anti-parallel to each other their emitted radiation inter- 
feres destructively. This implies that the system changes 
from superradiant behaviour around the primary reso- 
nance to a subradiant behaviour around the secondary 
resonance. Figure |4] displays the electric near-field and 
the dipole orientation at the two resonances. An in- 
tuitive explanation of the secondary resonance can be 
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FIG. 4: Electric field lines in the near field of two dipoles 
driven on the primary resonance (a), and the secondary res- 
onance (b) near a metal interface at z — 0. The bold arrows 
indicate the dipole location and orientation of the emitters. 



given as follows. For most driving field frequencies the 
coupling between the x and z components of the os- 
cillators is weak. However, at the resonance frequency 
of sub-radiant modes this coupling can transfer most of 
the energy of the oscillator from the z to the x compo- 
nents because the latter only lose very little energy energy 
through radiative emission and thus decay slowly. The 
secondary resonance therefore has to be narrow because 
its linewidth is inversely proportional to the large decay 
time of the x-oscillations. 

To improve our understanding we compare in Figs.lsFa) 
and (b), a numerical evaluation of the oscillation ampli- 
tudes for two emitters near the metal film with a cor- 
responding calculation for two emitters near a perfect 
mirror. The latter case can be described by taking the di- 
electric constant €2 of the metal film to be infinite, which 
corresponds to a perfect conductor. While the mirror im- 
ages explain the qualitative features very well, the nar- 
row resonance is greater in the case of the mirror. We 
attribute this difference to the existence of surface plas- 
mons, which are not present in the case of a perfect mir- 
ror. The dispersion and decay of surface plasmons cause 
the emitter energy to be dissipated quicker and this is 
represented by a corresponding decrease and broadening 
in the narrow resonance. 

With increasing distance from the interface the narrow 
resonance falls off exponentially because of the decreasing 
strength of the reflected near field at the interface. Pre- 
sumably for the same reason the resonance frequencies 
are shifted by a larger amount near the metal film. 



rate, respectively. Using the results of appendix [A| again 
we can transform Eq. ^ into a transcendental equation 
for the collective parameters uj and F, 



[{r^iujf^ujl]x^^\t) 
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G{r^^\r^'\uj-ir). (5) 



For optical emission it is safe to assume that F <C cj, 
i.e., the duration of the emitted light pulse is much longer 
than the inverse frequency. Because G{r^^\r^'^\uj)^ is a 
meromorphic function of Q in the lower half plane, and 
because its poles are determined by the properties of the 
metal interface and are unrelated to F, we can approxi- 
mate G(r, r^^\uj — iV) by G(r, r^^\uj — ze), where e > 
is infinitesimally small. This reduces the transcendental 
equation ([5| to a linear eigenvalue problem 
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(0) = 0, 



(6) 

with frequency shift S = ujo—uj <^ ujq. The decay param- 
eter F and the frequency shift 6 relate to the imaginary 
and real part of the eigenvalue problem (|6| , respectively. 
Generally, F characterizes the collective decay rate of the 
all oscillators. 



A. Decay rate 

It is instructive to first study the decay a single oscilla- 
tor x{t) at position r. In free space the Green's function 
takes the form [271 



G{r,r,uj) = [G' 



GttcJ 



(7) 



The real part G' is related to the Lamb shift of atomic 
resonance lines. It is formally divergent and would need 
to be renormalized [28]. However, because we are only 
interested in the decay rate, we can ignore the line shift 
and assume it is absorbed in the definition of the detun- 
ing. Because the Lamb shift is much smaller than the 
optical resonance frequency this will result in an excel- 
lent approximation for the decay rate. As in free space 
the matrix in Eq. m^ is proportional to the identity the 
corresponding eigenvalue problem is trivially solved and 
yields the well-known decay rate of a single oscillator. 



V. COLLECTIVE DECAY OF CLASSICAL 
EMITTERS 

We model the collective decay of initially excited os- 
cillators near the interface by assuming that there is 
no driving field and that the emitters perform a simple 
damped harmonic motion x^^\t) = ic^-^^(0)e~*'^^~^^, with 
UJ and F the collective oscillation frequency and decay 
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(8) 
(9) 



The decay rate near the interface can be derived in a 
similar way as long as F is much smaller than the opti- 
cal frequency and varies little for frequency variations on 
the order of the Lamb shift. These assumptions should 
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FIG. 5: The decay rate of a single oscillator placed within 
the vicinity of a thin metal film, at a separation distance of 
Az. the decay rates are given for a dipole oriented parallel 
(solid) and perpendicular (dotted) relative to the plane of the 
interface. 



be satisfied as long as the emitter is not too close to 
the interface. We remark, however, that these approx- 
imations are not universally valid and fail in the case 
of photonic band gap materials [29 , for instance. The 
Green's function G^j^{r^r^uo) evaluated at a single point 
r near the interface is diagonal in Cartesian coordinates 
(cf. App. [b|. We therefore can use Eq. ([8| to find the 
decay rate. 

Fig. [5] depicts the decay rate of an emitter as a func- 
tion of the distance Az of the emitter from the metal 
film (cf. Fig. fl]). It demonstrates that the decay rate is 
enhanced (suppressed) if the emitter oscillates perpendic- 
ular (parallel) to the interface, respectively. This can be 
understood using the concept of mirror images with each 
emitter considered as an electric dipole (see App.ffl. If 
a dipole is oriented perpendicular to the interface, its 
mirror image is in phase and thus can enhance the emis- 
sion of the dipole. On the other hand, dipoles that os- 
cillate in the plane of the interface have mirror images 
that are 180° out of phase so that the radiation reaction 
causes damping of the dipole. Hence, if the mirror images 
were composed of real charges these two situations would 
just correspond to superradiant and subradiant collective 
emission of two emitters, respectively. 

We now turn to the collective decay of A^ > 1 emitters. 
In this case there are 3A^ eigenstate solutions of Eq. ^. 
Fig. [6] displays four out of 24 numerically determined 
eigenmodes for A" = 8 which each are examples for a par- 
ticular collective behaviour. In free space a state where 
all emitters oscillate in phase and in the same direction 
is superradiant. The presence of the interface breaks the 
axial symmetry and oscillations in the plane of the in- 
terface are suppressed because the corresponding mirror 
images trigger subradiant behaviour. Fig. JGJ^a) shows the 
superradiant mode where the emitters oscillate in phase 
along the z-axis. In this case their mirror images are also 
in phase and enhance the radiation intensity. Fig. [6Fb) 
shows a subradiant mode where the emitted fields of the 
individual emitters interfere destructively. Fig. IgFc) and 
(d) display a very different state which we call suppressed 




FIG. 6: Four different collective eigenmodes: (a) superradi- 
ant mode, (b) subradiant mode, (c) and (d) are suppressed 
superradiant modes for which all real dipole moments are in 
phase but all mirror images are out of phase. Shown is the 
orientation of the dipole moment for real emitters (above the 
interface) and their mirror images (dashed, below the inter- 
face) . 



superradiant states. All oscillators are in phase and in 
free space it would be a superradiant state. However, the 
mirror images are out of phase so that the intensity of the 
emitted radiation is strongly reduced as compared to the 
superradiant state. In Fig. IsFa), which will be discussed 
below, we display the collective decay rate of these modes 
as a function of the number of oscillators. 



B. Far-field radiation 

A simple way to learn about the dynamics of a collec- 
tion of emitters is to observe their emission pattern in 
the far field. The radiation intensity is determined by 
the (time averaged) Poynting vector 



S{r,t) 
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(10) 



The electric field of each emitter is determined by Eq. ^ 
and the total field E{r^t) is the superposition of the in- 
dividual fields. To evaluate Eq. ([2| in the far field we 
have to compute the radiative Green's function (Bl) in 
the far- field limit, /cr ^ 00, which can be accomplished 
using the method of stationary phase described in Ap- 
pendix [C] 

For a source in the region z' > d and observation of 
the field at position r in the region 2: < 0, the dyadic 
Green's function can be approximated as 

27r \ r r J 



X G^jy{ki sin (j) COS 0^ ki sin (j) sin O^uj; z'\ 



(11) 



wi th th e unwieldy coefficients Gy,y defined in Eqs. (Bll) 
- (B19). A similar expression can be derived for an 




FIG. 7: The far- field emission pattern of 8 oscillators in the 
superradiant mode at a distance of l\z — 200nm from the 
interface (solid) and in free space (dashed). Shown is the 
ratio of the radiation intensity S and the maximum intensity 
S{) produced by a single emitter in free space. The emitters 
are placed parallel to the horizontal axis. 



observation point in region 3. The direction from 
the emitters to the observation point is given by 
(sin (\) cos ^, sin (\) sin ^, cos ^) . The radiation profiles are 
determined from the positive frequency components of 
the electric field in the far field. For emitters that are 
equidistant from the interface, the Poynting vector can 
be represented as 
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X \Gf^u{ki sin ^ cos ^, ki sin (j) sin ^ uj ; z )| 
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-ik sin 0(cos ^a^'-^-'+sin 9y^'^^) —Tt (i) 



(12) 



with j = 1,3 for the observation point in region 1 or 3, 
respectively. 

The last term is the relevant phasing term which de- 
termines the A/'^-gain in intensity that is associated with 
superradiance. The superradiant decay modes must be in 
phase, i.e., /csin0(cosfe^*^ +sin^?/^*^) + arg(x}i ) = ni27r 
for each emitter i, then the sum becomes proportional 
to N'^ . In other words, constructive interference between 
the emitted radiation of all emitters is one condition for 
observing superradiance. 



The transmitted radiation profile (12) for the super- 



radiant state in the plane of the emitters {9 = 0) is il- 
lustrated in Fig. [7| for a linear arrangement of 8 emit- 
ters periodically separated within a distance of A/2, with 
A = 800nm and u = 27rc/A. The dashed curve shows 
the emission pattern in free space. For perfect superra- 
diance its maximum intensity (along the horizontal axis) 
should scale like S = AN'^Sq, with ^o is the maximum 



intensity for a single emitter in free space and with the 
superradiance prefactor A = 1 in free space. Because the 
outermost oscillators are a distance A/2 apart so that 
their emitted fields are out of phase, the N'^ scaling is 
only roughly fulfilled. 

For emitters near the interface (solid curve in Fig. [tI 
surface plasmons lead to the narrow but extremely large 
peaks in the lower half of Fig. [7] at the plasmon reso- 
nance angle with a width of about 0.01 radian. Emission 
into the radiative modes around this peak is superra- 
diant (see below) with A ^ 240. The value of A de- 
pends on the separation distance of the emitters from 
the interface and the Fresnel coefficient T^2^ defined in 
Appendix [B] A will decrease exponentially as the sepa- 
ration distance between emitter and interface increases 
due to decreased near field coupling of the emitter to 
surface plasmon modes. For given dielectric constants 
Si and emission wavelength A there is an optimum film 
thickness for which A is maximized. The values used 
in this paper correspond to an optimum film thickness of 
56nm for A = 800nm. Hence, the radiation from multiple 
emitters transmitted into the surface plasmon resonance 
angle (the spikes in Fig. ItI) combines the enhancement A 
and the N'^ gain of superradiance. For the case N = 8 
the peaks have a maximum of almost 15000 times the 
maximum intensity So for a single emitter in free space. 



C. Suppressed superradiance 

We conclude this section with a discussion of the main 
indicators for superradiance in the superradiant and the 
suppressed superradiant states. In Fig. ^a) we display 
the scaling of the collective decay parameter with the 
number of oscillators, which is a measure for how strongly 
light emitted from one oscillator can drive emission from 
another oscillator. Because of the growing size of the 
arrangement of oscillators, which varies from 50 nm for 
A^ = 2 to 550 nm for N = 12, we expect superradiant 
phenomena to decrease with N. This is indeed the case 
for the superradiant state (circles), but for the suppressed 
superradiant states collective decay (F ~ N) (triangles 
and squares) is preserved. 

Interestingly, we observe the opposite situation with 
respect to the N"^ dependence of the peak intensity at 
the surface plasmon resonance angle, which is shown in 
Fig. [SJ^b) . For the superradiant state the peak intensity 
can well be described by 5 ~ AN'^Sq^ indicating con- 
structive interference of the emitted radiation despite the 
growing size of the sample. For the suppressed superradi- 
ant states the ability to constructively interfere depends 
on the dipole orientation: the N'^ scaling is significantly 
affected only if the dipoles are aligned as in Fig. p[c). 

Clearly the two suppressed superradiant states fulfill 
the criteria for superradiance as well, or better, than the 
superradiant state itself. We use the notion "suppressed" 
to characterize their behaviour because their decay rate 
is reduced by a factor of 6 and their peak intensity by 
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this direction phase variations have a very pronounced 
effect. For the suppressed superradiant state [6jd), as for 
the superradiant state, most radiation is emitted into a 
plasmon that propagates in the ^/-direction, i.e., paral- 
lel to the alignment of emitters. Therefore their emitted 
radiation fields can still constructively interfere. 



FIG. 8: (a) Scaling of the collective decay rate with the num- 
ber N of oscillators. Circles correspond to the superradiant 
mode, squares and triangles to the suppressed superradiant 
modes displayed in Fig. Ig] (c) and (d), respectively, and dia- 
monds to the subradiant mode displayed in Fig. p] (b). The 
distance between adjacent oscillators is 50 nm. Each oscilla- 
tor is a distance Az — lOOnm away from the interface. In 
(b) the scaling of the maximum intensity of the peak at the 
surface plasmon resonance angle with the number N of oscil- 
lators. In the superradiant case the results are reduced by a 
factor of 10 for presentational purposes. 



a factor of 25 as compared to the superradiant state. 
We can understand suppressed superradiant behavior by 
considering each dipole and its mirror image as one en- 
tity. The dipole moment of this entity vanishes, but it 
can emit (a weaker) quadrupole radiation. Because all 
N quadrupoles are in phase their emitted fields inter- 
fere constructively and the intensity scales like N'^. Fur- 
thermore, similar to how a neighbouring dipole in phase 
can increase the decay rate, a quadrupole can drive the 
neighbour when in phase. Hence, the suppressed super- 
radiant state can be considered a superradiant state for 
quadrupoles. 

The different behaviour of the states of Fig. [6] can 
be understood by considering the mirror images as real 
charges. For suppressed superradiant states the spatial 
size of each quadrupole, 2Az = 200nm, is larger than 
the distance Ax = 50nm between adjacent quadrupoles. 
Hence, they interact with their respective near field so 
that cooperative decay can still grow like N. On the 
other hand, for the superradiant state each dipole inter- 
acts with the far field of other dipoles so that the devia- 
tion from the linear scaling with N is more pronounced, 
albeit the overall interaction energy is much larger. 

The 7V^ scaling of the peak intensity can be under- 
stood by taking into account that the maximum peak 
appears at the surface plasmon resonance angle, i.e., in a 
direction that is different from the axis along which the 
oscillators are aligned. Generally, surface plasmons can 
only be generated in the direction of the polarization of 
the emitted radiation because it is the latter's electric 
field that generates electron density oscillations. For the 
suppressed superradiant state p[c) there are significant 
deviations from the N'^ law because the dipoles oscil- 
late in the plane of the interface parallel to the x-axis. 
This means that most radiation is emitted into a surface 
plasmon mode that propagates in the x-direction, but in 



VI. CONCLUSION 

In this paper we give a detailed account of how super- 
and sub-radiant emission by oscillating dipoles is affected 
by surface plasmons. It was shown that the interface gen- 
erates an indirect coupling between the emitters through 
light that is refiected by the interface. In the superradi- 
ant mode the peak intensity of the emitted radiation is 
two orders of magnitude larger than a similar arrange- 
ment of emitters in free-space. The decay rate of the 
system is increased by the presence of mirror dipoles at 
the interface. Applying a driving field to the emitters 
can induce super- and sub-radiant emitter modes by a 
suitable choice of frequency. 

The collective decay eigenmodes of initially excited os- 
cillators showed that certain modes, which naively would 
be expected to behave in a sub-radiant fashion, are ac- 
tually superradiant modes, albeit with a strongly sup- 
pressed overall intensity. These suppressed superradiant 
modes are also a consequence of the additional coupling 
between oscillators that is generated by the interface. An 
intuitive picture explains this effect as superradiance of 
quadrupole radiation. 

Our results will help to improve the characterization of 
the radiative response of entangled atoms or ions in quan- 
tum information and of surface-enhanced spectroscopy of 
molecular excitations. A detailed study of photon entan- 
glement in the presence of surface plasmons is currently 
under investigation. 
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Appendix A: Derivation of the dynamical equation 

The propagation of an electric field E{r^t) in the pres- 
ence of multiple interfaces is solved using a Green's Func- 
tion approach. The electric field due to a current source 
of j(r^t) is determined by the well known differential 



wave equation 



-sir) (^ 



F V X Vx J E{r,t) = -no^^j{r,t), 

(Al) 
where €{r) is the complex relative permittivity of a di- 
electric. We assume here for simplicity that the permit- 
tivity of the interfaces does not depend on the radiation 
frequency, which is the case for the interface that we are 
studying [35 . The constants c and /io are the speed of 
light and permeability of free space respectively. The so- 
lution to Eq. ( Al ) in terms of the dyadic Green's function 



G(r,t]r',t') is 



E{r,t) = -fioJJdt'd'r'G{r,t;r\t')-^j{r\t'), (A2) 

We assume that each emitter consists of a charge q (the 
"electron" ) that performs small oscillations of amplitude 
x^^\t) around a fixed position r^^^ at which a charge —q 
(the "nucleus") is placed. The current density j{r\t') 
produced by emitter j can then be written as 



j{r\t') = qS{r' -r^^^)x^^\t') 



(A3) 



This current density contains the complete information 
about the radiation field produced by both the electron 
and the resting nucleus and determines the radiative elec- 
tric field through Eq. (A2). In addition one generally 



needs to include the electric Coulomb field of a charge dis- 
tribution. However, because the oscillation amplitudes of 
the electron are small and the total charge of the emit- 
ter is zero, each oscillator essentially describes a point 
dipole with dipole moment p = qx{t). We then can ig- 
nore the Coulomb contribution so that each emitter can 
be thought of as a n os cillating point dipole. Inserting 
Eq. (|A3| into Eq. ([A2| then yields the electric field ^ 
produced by emitter j. 

The radiative Green's function G{r^t]r' ^t') is a re- 
tarded solution to 



-sir) (^ 



Vx Vx G{r,t;r\t') 
= IS{r-r')5{t-t'). (A4) 



The Green's functions in frequency domain and Cartesian 
coordinates are derived from the equation 



d d 



S^^ ( A 



-e(r,cj) 



G^^(r,r',cj) 
S^J{r-r'). (A5) 



where the indices indicate the appropriate cartesian co- 
ordinates. Considering that the sources are located in 
region 3, the solution to Eq. (A5) is [24|, 



G^y{r,r',uj) 



2 d d 



ik^\r-r'\ 



AttW ■ 



TTr (A6) 



where k^, = >/^^- "^^^ Green's functions corresponding 
to the refiected and transmitted fields are obtained by 
solving the homogeneous equation 



d d 



J^^ ( A + es ^2 



G^,(r,r^a;) = 0. (A7) 



The solution to these equations with applied boundary 
conditions are found in Appendix [B] 

Let us now assume that emitter j oscillates according 
to 



x^^\t) = x^^\0)e- 



-iujt-rt 



(A8) 



for some frequency oo and a real, positive decay parameter 
r <C cj. Eq. Q then becomes 

/CX) 
dt'G{r,r^^\t-t') 
-CX) 

x(r + iw)2a;(^')(0)e-*'"*'-r*' 

xG{r,r^^\oj-ir)x^^\0) , (A9) 

where G{r,r^^\oj — IT) is the analytic continuation of 
the temporal Fourier transform of G{r,r^^\t) into the 
lower half plane. 

Noting the pole at wi = to — iT and applying calculus 
of residues, the electric field is then 

^(i)(r.,i) = -iiQqG{r,r'^^\LO-iT){T+iu)fx'^^\Q)e-'^*-'^K 

(AlO) 
Inserting Eq. ( |A10| ) and Eq. ( [AS] ) into Eq. Q gives 



ml 



m 

N 



X "^G{r^^\r^'\uj-ir) 

i=l 

X {r^iujfx^'\t)\ (All) 



Appendix B: Green's Functions of a 3-layer 
dielectric 



Using the method of Ref. [30^ (See also Ref. [3T1) the 
dyadic Green's functions of a multiple layered planar di- 
electric are solved. Due to the translational invariance of 
the problem of a planar dielectric interface, it is useful 
to decompose the Green's function into transverse and 
normal components through the Fourier transform 

G^,{r, r', uj) = j ^e^K-(r,-r'^)G^,{kp, to; z, z'), 

(Bl) 
where kp is the wave vector component tangential to the 
interface and Vp is the corresponding position component. 



For sources located in region 3 {z' > d) and z > d the 
elements of G^j,{kp^ uo]z^z') are 

,{K,z')) (B2) 



'^^y^3,2^yy\ 



^xy 



fh'T 1^' 



_ is^xi^y jp^(z-d) ( 5TM 



L.2 



( -^3,2 ^xxyf^p^ Z ) 



+ il3 2 Z/yy[Kp^ Z 



^<J^xz — , -^3,2 "^xzyl^p^ ^ J^ 



^yx — ^x 



^yy — ^xxi l^x "^"^ l^y 
^yz ^xzi f^x "^ f^y 



G. 



f^X 



-^3,2 "^xxyf^p^ ^ j6 



'zy n 



Ps 
and for z < the terms are 

^xx = 1^^ l^a?^3,l "^xxyf^p^ ^ ) 



+k. 



k h 



^xy — h,2 ^ \ 3)1 "^XXX^^pT ^ ) 






^yx ^xy 

^yy ^ ^xx") f^x "^ i^y 

^XZl 

I^X Q 

roil p; 



^yz — ^xz-} i^x "^ i^y 



Ps 



Ps 



where 



^zz — -o ^3 1 ^xz\i^p-) Z j6 

Ps 



ips 



ZiyyyrZp^ Z ) — _ ^ 6 



Zrv'T) 



(B3) 



xz — ,^ -^^3.2 ^xzy^pi "^ y- v-'^'^/ 



(B5) 
(B6) 
(B7) 

(B8) 



G.„ = ■fR™Z,,{k„z')e'^^^^-''^ (B9) 

53 



(BIO) 



■^ySjPyy{kp,z')) (Bll) 



Sj^Zyy{k„z')) (B12) 

53Vr^..(fcp,^'e-*''^^) (B13) 



(B14) 
(B15) 
(B16) 



zx — r^ ^3,1 ^a^xi^p?^ j6 (l^-'-'j 



G,„ = ^5™Z,,(fep,^')e""'^" (B18) 



(B19) 



^..(fep,^') = S^'"'""''' (B20) 

'^ = J_il3s\d-z'\ 
2ft 



(B21) 



^..(fep,^') = ge^'^^M-^'l, (B22) 



with /3^ = Jkf — k'^ and /c^ = y^^- 

The transmission and reflection of radiation at the in- 
terface is characterized by the generalized Fresnel reflec- 
tion and transmission coefficients. Consider a series of 
interfaces, where i increases for z increasing in 2: > 0. 
For reflection in region i, from an interface between re- 
gions i and i — 1, the generalized reflection coefficient 
is 



R 



,TM/TE 



^TM/TE 5TM/TE .2/3,_i(d,_ 
^i,i-l '^ ^i-l,i-2^ 

^ ^TM/TE5TM/TE ,2/3._i(d.. 



The generalized transmission coefficients are 



nTM/TE 



n 



rpTM/TE if3^^drr 

^m,m-l ^ 



-di-2) 

(B23) 

,-1) 



m=j-\-l 



^TM/TE^TM/TE i2/3„_i(d„_i-d„_2) ' 



1 



(B24) 
These are expressed in the usual Fresnel coefficients 



R 



TM 
iA—l 



pTE 



nTM 



^TE 

■^iA-1 



^i- 


-ift-eift- 


-1 


ft 


_lft+£ift_ 

-ft-i 


-1 


ft + ft-i 

2£ift-i 




£i- 


-ift+£ift- 
2ft 


-1 



P^ + P^- 



(B25) 
(B26) 
(B27) 
(B28) 



Appendix C: Stationary phase method 



To evaluate the two-dimensional integral in Eq. ( |B1 ) 
we use the stationary phase method (see, e.g., Ref. [32 J. 
There are two contributions to the phase of the integrand: 
the Fourier transform term and the phase of the Green's 
function itself, which depends on the phase of the Fres- 
nel coefficients and a phase factor that depends on the 
position z' of the source. Because the film thickness d is 
much smaller than A the Fresnel coefficients vary slowly 
with kp. For this reason we need only to consider the 
position of the source. For brevity we here will only dis- 
cuss the Green's function for the observation point r in 
region 1 and the source position r^ in region 3. Because 
of the general form of the Green's function in this case 
(see Eqs. (Bll) - (B19)) integral (Bl) can generally be 
written as 



G^^(r,r',cj) = / 



Qj l\"n 



(27r)2 



^ikp-{r-r)^-i(3iz^i(33\d-z'\ 



(Cl) 



where the coefRcients D ^^{kx , ky , to) are related to the 
Fresnel coefficients and can be formally defined as 
D^,{kx,ky,uj) = G^,{kp,uj;z,z')e'^'^e-'Md-z'\^ ^he 

phase term of interest is then 



^ = kp-{r-r')-piz + l33\d-z'\. 



(C2) 
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Recognizing that z' > d, and that kp = (kx^ky^O) is the 
same in region 1 and 3 we get 

^ = kx{x - xO + ky{y - y') - ^z + p^{z' - d) (C3) 

The phase is then expanded in powers of kir. To do so 
we introduce the quantities 



^x — 



y 



r 

x[ 

r 



ki ki 



z 



,2x1/2 



^z ~ ~ ~ ^ ^ ^x ~ ^y ) 

/'l (-1 2 2\l/2 

m = —-;—= —{I — p — q ) ' 
fci 



(C4) 
(C5) 
(C6) 
(C7) 
(C8) 
(C9) 



fci 



2 \ 1/2 

I -f^-^M . (CIO) 



where {sx^ Sy^Sz) determine the direction of the observa- 
tion point relative to the source, and p, q are re-scaled 
integration variables. We can then express the phase as 

<|) = kir {p{sx - s'^) + q{sy - Sy) + msz + rn's^) . 

(Cll) 
The integral is evaluated in the limit of large kir. The 
stationary points are determined by setting the deriva- 
tives d^/dp and d^/dq to zero which yields 



^^-^^ -P- and -'-''^ - ^ 



Sz + s'^ m 

z z rn' 



. (C12) 



Because r ^ r' we generally have s[ <^ Si. This is not 
true if we observe the far field very close to the z-axis (so 
that Sx ^ s^^) or very close to the plane of the interface 



(so that Sz ~ s'^). Ignoring these special cases we can 
simplify the stationary points to 



^ = P and ^ = ^, 

Sz m Sz m 



(C13) 



which is the usual result of the far field approximation 
that k and r are parallel. 

To determine the contribution of the critical points to 
the integral the second order partial derivatives at the 
stationary points must be evaluated. These evaluated 
derivatives are 



dx'^ 

dy'^ 
dxdy 




~x'^y 



(C14) 

(C15) 
(C16) 



If we parametrize the stationary point as p = I^l = 

sin cos ^, and q = -^ = sine/) sin ^, the Green's func- 
tion in stationary phase approximation becomes 



Gu 



Ar,r' ,UJ) = -^ ( — -] ^ikir ^-ikisincl){cosex' +siney') 

27T \r r J 

^ ^-i(3iz^i(33\d-z'\ 

X Dj^iy {ki sin (j) cos 6>, ki sin (j) sin 6>, u). 

(C17) 



The terms I3j = kz^i are evaluated at kp = ki sin0, such 
that Pj - 



k'^ — kf sin^ ^. The effect of the stationary 

phase method on the integral is a leading term that is 
proportional to z. This results in the vanishing of terms 
near the boundary with the exception of those tied to TE 
radiation. Expressing D^j^ in terms of G^jy then yields 
Eq. ([TTl). 
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